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ABSTRACT 

HM Cancri is a candidate ultracompact binary white dwarf with an apparent orbital period of only 
5.4 minutes, as suggested by X-ray and optical light-curve modulations on that period, and by the 
absence of longer-period variability. In this Letter we present Keck-I spectroscopy which shows clear 
modulation of the helium emission lines in both radial velocity and amplitude on the 5.4-minute period 
and no other. The data strongly suggest that the binary is emitting Hei 4471 from the irradiated face 
of the cooler, less massive star, and Hell 4686 from a ring around the more massive star. From their 
relative radial velocities, we measure a mass ratio q — 0.50 ±0.13. We conclude that the observed 5.4- 
minute period almost certainly represents the orbital period of an interacting binary white dwarf. We 
thus confirm that HM Cnc is the shortest-period binary star known: a unique test for stellar evolution 
theory, and one of the strongest known sources of gravitational waves for the Laser Interferometer 
Space Antenna (LISA). 

Subject headings: binaries: close — X-rays: binaries — gravitational waves — stars: individual (HM 
Cancri, V407 Vulpeculae) 



1. INTRODUCTION 



Two interacting binary stars have been discovered 
which appear to have orbital periods shorter th an ten 
minut es: V407 Vul, with a period of 9.5 minutes (|Motch 

of 



et al. 1996 Cropper etal. 



1998|), and H M Cnc, with a 



period of 5.4 minutes ( jlsrael et al. 2002). The uniquely 
short period of 5.4 minutes, if it is the orbital period, 
implies that HM Cnc must have formed from two white 
dwarfs, driven together as a result of gravitational-wave 
radiation. It may currently be experiencing stable mass 
transfer through Roche-lobe overflow. Because it is po- 
tentially so extreme and unique, substantial effort has 
been put into unveiling HM Cnc's true nature, but as 
yet without conclusive results. 

The key observational data show that: (1) there is no 
evidence for variability on periods other than 5.4 a nd 9.5 



minutes in HM Cnc a nd V407 Vul, respectively (Ram 
say et aL]|2000 2002a); (2) the optical flux maxima lead 



(Barros et al. 


200' 


and (3) the periods are decreas- 


ing in both systems (Strohmayer 2002, 2003 


2004 


2005 


Hakala et al.||2003 


2004||lsrael et al.||2004). 



Vul and HM Cnc. One of them, the 'Intermediate Polar 
(IP)' model, predicts that these systems are not in fact 



ultracompact binaries but have rather mundane orbital 
periods of several hours. The ultrashort pe riods then 
represen t the spins of magnetic white dwarfs ( Norton et 
aL||2004 1. The X-rays as well as the variable optical flux 



originate from the accretion flow crashing onto the mag- 
netic poles of the magnetic white dwarf during part of the 
accretor's spin cycle. The spin-up of the X-ray and op- 
tical periods may be expected since the magnetic white 
dwarf is accreting matter of high specific angular mo- 
mentum. The absence of variability on the much longer 
orbital period could be explained if the orbital planes 
of both systems are viewed exactly face-on. The ob- 
served phase offsets between th e optical and X-ray light- 
curves are difficult to explain (Norton et al. 20041, and 
the emission lines are unusually weak for IPs; in V407 



Vul there appear to be no emiss ion lines at all (Ramsay 



Vul t nerc appear t o pe no emissi 
et al.||2002b[ Steeghs et al.||2006[ ) 
The second, 'Unipolar Inducto 



polar Inductor (UI)' model, is essen- 
tially a more energetic version of the Jupit er-Io system 
( |Wu et al.|20"02{|DairOsso et al.|2006[|2007| ). A magnetic 
white dwarf is orbited by another, non-magnetic one; the 
magnetic field induces an electrical potential across the 
non-magnetic white dwarf which sets up currents along 
magnetic flux tubes connecting the two stars. Ohmic 
dissipation of these currents in the flux tubes' footpoints 
on the magnetic star gives rise to the X-rays. In this 
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model the 5.4-minute period is the orbital period, but 
the two stars are detached. Since the binary loses an- 
gular momentum due to gravitational-wave radiation, it 
is expected to evolve towards shorter orbital periods, as 
observed. The main problem is the offset between the op- 
tical and X-ray flux, which requires the footpoints on the 
magnetic star's surface to be almost 90 degrees a head of 
the orbiting non- magnetic star in both systems (Barros 
et~aL|2005| 2007| . 

I'he third, 'AM CVn' model, has a Roche-lobe filling 
white dwarf losing mass to a more massive w h ite dwarf 
dMarsh fc Steeghs|2002||Ramsay et al.|2002b| or |Cropper ' 
|et al.||l998| tor a magnetic version), i'he orbit al period 
is 5.4 minutes in the case of HM Cnc, and the accretion 
stream hits the accretor directly without forming a disk. 
The phase offset between optical and X-ray flux is natu- 
rally accounted for, since the accretion stream will deflect 
and hit the accretor off-axis. The absence of longer pe- 
riods is also expected, but the observed decrease of the 
periods in HM Cnc a nd V407 Vul is considere d problem- 
atic for this model (Strohmayer 2002 2004), al though 
solutions to this problem have been put forward (Deloye 
|fc Taam|200~6l|D'Antona et al.|2006| Roelofs & Deioye in 
preparation). 

This Letter describes our search for kinematic evi- 
dence, using phase-resolved spectroscopy of HM Cnc, for 
or against these different models. Section §[2] describes 
our observations and data reduction. We present our 
results in §[3] and discuss their implications in §|4] 

2. OBSERVATIONS AND DATA REDUCTION 

After unsuccessful attempts in 2005, 2006 & 2007 to 
obtain phase-resolved spectroscopy of HM Cnc, we finally 
met acceptable weather conditions at Keck-I on January 
25 & 26 and March 30 of 2009. We took spectra using a 
1.5" slit and the 600 grooves/mm grism on the blue side 
of LRIS, th e Low-Resolution Imaging Spectrograph ( Oke 
et al.|[l995 ), for an effective spectral resolution of about 



300 km/s. In order to resolve the 321.5-second period in 
HM Cnc, we limited the exposure times to 60 seconds. 
The blue CCD, with a plate scale of 0.135"/pixel, was 
binned 4x4 pixels to reduce the read-out noise while 
still sampling resolution elements with approximately 2 
(binned) pixels. The binning furthermore reduced the 
dead-time between exposures to 27 seconds. 

Data reduction was done using standard iraf tasks. 
To optimally extract the very faint individual spec tra 
(HM Cnc has op tical magnitudes U » 19.6, B « 20.7;[lT" 
rael et aL]|2002[) , reference aperture profiles were created" 
from averaged blocks of typically ^30 successive frames 
to increase the photon statistics and remove cosmic rays 
from the profile. Spectra were grouped together based on 
the requirement that spatial center and profile (seeing) 
on the CCD be (near-)identical. Wavelength calibration 
was achieved using HgCdZnNeAr arc exposures. A fit of 
30 arc lines left 0.35A root-mean-square residuals. No 
detectable arc drift occurred during the night, and the 
drift between the two nights was well below one pixel. We 
thus applied a common dispersion solution to all spectra, 
only applying a small shift to the first night's spectra to 
tie the strong [O i] 5577 sky emission lines together. The 
relative wavelength calibration of all spectra is estimated 
to be better than 0.1 A. 

The time stamp and velocity scale of each spectrum 



were transformed to the solar system's barycenter. Time 
stamps were checked several times during our run using 
Universal Time clocks available online and estimated to 
be correct to one second (or better). The spectrophoto- 
metric standard star Feige 34 was used for calibrating the 
instrumental response; variable transparency and seeing 
precluded an absolute flux calibration. Table [T] summa- 
rizes our observing log. 

3. RESULTS 
3.1. Average Spectrum 
The average Keck spectrum of HM Cnc is shown 



in Fig. [I] As previously reported (Israel et al. 2002[ ) 
the spectrum is dominated by ionized helium emission 
lines. Their full-width at half-maximum (FWHM) of 
^2500 km s -1 is well resolved by the 300 km s -1 reso- 
lution of our spectra. Lines of neutral helium are also 
present but much weaker. Our higher-quality spectrum 
confirms earlier findings that the even-term transitions 
of the Hell Pickering series are stronger than t he odd- 
term ones (s ee Fi g. [T|, which was interp reted by |Norton| 



et al. (2004) and Heinsch et al. (2007) as evidence tor 



the presence of hydrogen. In particular the absence of 
Hen 4200 in our spectrum is striking, and most easily 
explained if one assumes that there is indeed hydrogen 
in HM Cnc (for simplicity, we will refer to He II Pickering 
lines in the remainder of this paper). 
In addition to the Nm 4640 line in the Bowen blend 



identified by |Israel et al.| ( |2002[ ), the strong N in 4379 line 
appears to be present to the red of Hell 4338. Apparent 
emission lines near 4600 A and 4800 A could be due to 
N iv, while a faint feature to the red of He I 4471 may be 
N in 4514. The presence of several nitrogen lines appears 
to be at odds with the suggestion (based on fits to X-ray 
spectra) that HM Cnc may have an unusual chemical 
abundance pattern and be un derabundant in nitrogen in 
particular ( |Strohmayer 2008). We note that nitrogen- 
rich matter would be expected due to the CNO-cycle if 
HM Cnc is an ultracompact binary, unless th e material 
experienced helium burning (| Yungelson| [2008 ) , and that 
nitrogen i ndeed appears to be abundant in most AM CVn 
stars (see Roelofs et al.||2009 ). 



al. 



3.2. Ephemeris and Period Derivative 

When we phase- f old our spectra using the ephemeris of 
( 2007 ), we observe continuum flux variations 
ltly (~30 degrees in phase) behind Barros et 
(|2007|'s photometry. This is not unexpected given 



Barros et al. 
that lag slig 



m particular , the uncertainty in the period derivative 



measured by Barros et al. (2007), which translates into 
a phase uncertainty of about this magnitude. 

By matching the phases of the variable <?-band flux in 
our time-resq l yed sp ectra with the g-band photometry of 
Barros et al. (2007), we can obtain a refined ephemeris 
for HM Cnc. This will not be a very precise refinement 
since our spectra only marginally sample the orbital pe- 
riod, and since atmospheric transparency and seeing were 
variable. However, because we have hundreds of spectra 
and the folding period is much shorter than the typical 
timescales for sky quality variations, we can still resolve 
this shift. The easiest way to correct the phase shift is 
to reduce the frequency v and its time derivative from 
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TABLE 1 

Log of our Keck/LRIS observations of HM Cnc. 



Date 


UT 




Spectra 


Weather 


Total S/N 








(60s) 




(at 4200 A) 


2009/01/26 


08:10- 


-12:30 


136 


1.1", absorption early 


78 


2009/01/27 


07:24- 


-13:14 


223 


1.1", clear 


115 


2009/03/31 


05:52- 


-06:50 


41 


1.3", clear 


41 




4000 4500 5000 

Wavelength (A) 

Fig. 1. — Average Keck-I/LRIS spectrum of HM Cnc, normalized 
to the continuum. The most obvious spectral lines are labeled. 
The Hell 3888 and 4026 lines may contain a contribution from 
coincident He I lines. 



Barros et al.| ( 2007| ) by approximately 1.3 a: 

to (barycentric TDB) = 53009.889943753 
z/(Hz) = 0.00311013811(10) 
dv/dt(Uz s- 1 ) = 3.57(2) x 10~ 16 



(1) 



where numbers in parentheses are the uncertainties in the 
corresponding number of last decimals. The new error 
on the frequency derivative is estimated as the error that 
gives the same phase uncertainty (at our epoch) as the 
error on the frequency. The error on the frequency has 
been kept the same. This ephemeris is used throughout 
this paper. 

3.3. Time-Resolved Spectrum 

The left panels of Fig. [2] show the time-resolved 
(trailed) spectrum of HM Cnc around the He I 4471 line, 
folded on the above ephemeris (eq. (fll)). The continuum 



flux is seen to fo llow broadly the same pattern as in Bar- 
ros et al.] ( |2007[ ) after aligning the phases: during a cycle 
the continuum flux increases slowly and decreases more 
rapidly. The scatter in our continuum fluxes, shown in 
the left-most panel of Fig. [2] is due to differences in the 
average sky quality between phase bins rather than due 
to noise in individual spectra. 

Clear spectral line variations are seen in the He I 4471 
line, in the form of an 'S-wave' feature that is visible for 
about half the period while it moves from red-shifted 
to blue-shifted wavelengths. The intensity of the S- 
wave in the Hei 4471 line follows the intensity of the 
variable continuum flux, suggesting that they originate 



from the same region. In the remaining panels of Fig. 
[2] the intrinsically variable continuum flux has been fit- 
ted and subtracted to highlight the line variability. Here 
the Hen 4686 line is plotted as well as the sum of the 
three strongest He 11 Pickering series lines (4100, 4338 & 
4859 A; see Fig. 111). The He 11 lines behave quite dif- 
ferently from theHei 4471 line. The Hen 4686 line 
is double-peaked and appears to wobble in anti-phase 
with respect to the Hei 4471 line. The He 11 Picker- 
ing lines have a strongly modulated, narrow component 
which may be moving from blue-shifted to red-shifted 
in phase with (and at the same radial velocity of) the 
Hell 4686 line. In addition the He 11 Pickering lines have 
a fairly constant, broad component with a FWHM of 
-2500 km s" 1 , which matches the FWHM of the (double- 
peaked) He 11 4686 profile. 

We me asure, from a linear bac k-projection Doppler to- 
mogram (Marsh & Home 1988), a radial velocity semi- 
amplitude of 390 ± 40kms _i for the S-wave feature in 
the Hei 4471 line (corrected for a —5% bias due to fi- 
nite exposure times, which we modeled using synthetic 
data) . The error was estimated from a large ensemble of 
Doppler tomograms calculated from bootstrap samples 
of our 400 spectra. For the double-peaked He 11 4686 line 
moving in anti-phase, we similarly measure a radial ve- 
locity semi-amplitude of 260±40kms _1 . The two peaks 
in the profile are red- and blue-shifted by 650±50 kms -1 . 

Finally, the spectra were phase-folded on trial periods 
of up to 6 hours, but no radial- velocity or intensity vari- 
ations in the spectral lines were seen on periods other 
than 5.4 minutes. 

4. DISCUSSION AND CONCLUSIONS 
4.1. The Nature of HM Cnc 

The observed radial- velocity variations in the spectral 
lines strongly suggest that HM Cnc has a 5.4-minute or- 
bital period. We favor the semi-detached binary white 
dwarf (or AM CVn) model, shown in Figure [3j It is the 
only model that can account for the spectroscopy pre- 
sented in this Letter. The IP model predicts an orbital 
period of hours rather than minutes, and the observed 
line ki nematics do not match predictions by |Norton et al.| 
(2004) for spectral line variability on the accretor's spin 
period (if any). We see the broad, fairly constant (and 
in the case of He 11 4686, double-peaked) He II emission 
as a clear signature of accretion, which is not predicted 
by the UI model. 

We conclude that HM Cnc is by far the shortest-period 
binary star known. Its direct progenitor was probably a 
detached binary white dwarf, since evolutionary models 
with initially non-degenerate donor stars cannot reach 
such short orbita l periods (e.g. Yungelson||2008 van der 
Sluys et al.|[2005 ). 

The observed decrease of HM Cnc's 5.4-minute period, 
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Fig. 2. — Trailed spectra of HM Cnc, folded on the ephemeris of Eq. uXJ) , and resampled into 11 phase bins. The left-most panel shows 
the continuum flux variations near the He I 4471 line (dots with error bars} as well as the line emission intensities, measured between ±800 
km/s relative to the continuum (triangles; scaled up by a factor five). In each but the left-most trailed spectrum the varying continuum 
flux has been subtracted to highlight the emission line variability. The middle trail shows the double-peaked Hen 4686 line moving in 
anti-phase with the He I 4471 line. The right two panels show the sum of the three strongest Hell Pickering series lines (4100, 4338 & 
4859 A; see Fig. [IJ, where in the right-most panel the average spectrum has been subtracted to highlight the time- varying component. The 
grey-scales are linear between the dimmest and brightest pixels, where the brightest pixels have the highest flux densities. 
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Fig. 3. — 'AM CVn' model of HM Cnc, with the donor star on the 
left, transferring mass to the accretor. The impact of the accretion 
stream into the accretor causes the X-rays. The irradiated face of 
the donor star is the source of the He I 4471 emission line, and the 
main source of the modulated optical emission. An equatorial belt 
or disk around the accretor is the source of the Hell emission lines. 
Masses M 2 = 0.27 M , Mi = 0.55 M are assumed (see text); the 
center-of-mass is at the origin. The orbital phase (from Eq. |ljl) 
would be ~0.35 as viewed from the bottom of this page. 



considered the main obstacle for the AM CVn model, 
may be expected naturally if it is a semi-detached binary 
white dwarf with this orbital period, due to the relatively 
long mass-transfer turn-o n timescales predicted for fairly 
low-mass helium donors ( Deloye fc Taam||2006[ Roelofs 
& Deloye in preparation). 



4.2. System Parameters 

Based on the radial velocities of the He I 4471 and He n 
4686 emission lines moving in anti-phase, we can estimate 
the ratio of donor to accretor mass q = M2/M1 in HM 
Cnc. The He I 4471 line seems to originate from the 
irradiated face rather than the exact mass center of the 
donor star, and so we must estimate a 'K-correction' (e.g. 



Mufioz-Darias et al. 2005) to obtain the donor's actual 
projected orbital velocity from the observed emission-line 
velocity. Without a K-correction we get an upper limit 
of q < 0.67±0.12. Assuming that the Hei line originates 
from the inner Lagrange point (the maximum conceivable 
K-correction) gives a lower limit q > 0.33 ± 0.06. The 




-1000 



1 000 



K x (km s ) 

Fig. 4. — Doppler tomograms of He I 4471 (gray-scale) and Hell 
4686 (contours). The (assumed) irradiation-induced He I 4471 
emission from the secondary star has been aligned with the pos- 
itive ify-axis. The gray-scale indicates the flux overdensity as 
a percentage of the continuum. Hell 4686 is emitted in a ring 
centered on the expected location of the primary star, for mass 
ratios q = 0.50 ± 0.13. Solid and dashed contours indicate where 
the He II 4686 flux rises above and drops below the average in the 
map, respectively. A bright spot occurs on the left side of the ring, 
approximately 110 degrees ahead of the secondary star. 



true value will be in between; assuming a flat probability 
density distribution between our upper and lower limits 
gives q = 0.50 it 0.13. 

Figure [4] shows a combined Doppler tomogram of the 
He I 4471 and Hen 4686 lines, showing the (putative) ir- 
radiated donor star and ring-like emission centered on the 
(proposed) accreting star. The mass ratio derived above 
can easily be estimated from this figure. A further inter- 
esting feature is the 'bright spot' in the He 11 4686 emis- 
sion, occurring in the lower left quadrant of the Doppler 
tomogram. The velocity vectors from the center-of-mass 
to the donor star and to this He 11 4686 bright spot make 
a 110 degree angle approximately. Kinematically, this 
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bright spot matches with a spot approximately on the 
side of the accreting star as seen from the donor, if the 
material in the bright spot is moving at roughly half the 
Keplerian (break-up) velocity near the surface of the ac- 
cretor (but depending on its mass) . An accretion stream 
impact sp ot at this location was proposed by |Barros et] 
al. ( 2007 ) to explain the relative phases of the optical and 
X^ray light-curves. Assuming M2 = 0.13 Mq, which is 
the minimum donor star mass corresponding to a Roche- 
lobe-filling degenerate helium object, an accretor mass 
Mi w 0.55 Mq is required to get an accretion stream im- 
pa ct spot sufficiently on the side of the accretor, as shown 
by Barros et al. (2007). Our mass ratio q = 0.50 ± 0.13 
implies that both the donor and the accretor have to be 
more massive than this to still produce an impact spot 
sufficiently on the side (but see Wood|20"09 for a possible 
relaxation of this constraint). 

The measured rate of change of the orbital period in 
HM Cnc (Eq. [T| also suggests that the donor star has 
to be significantly more massive than its fully-degenerate 
value, if one assumes that the rate of change matches the 
secular rate set by gravitational-wave radiation (neglect- 
ing mass transfer). Given q ss 0.50, masses M 2 = 0.27, 
Mi = 0.55 would be required to yield the observed 
period derivative. We note that these are quite typ- 
ical masses for binary white dwarfs a ccording to the- 
oretical models (Nelemans et al. 2001). Assuming in- 
stead M2 = 0.13 Mq, which is again the minimum mass 
of a Roche-lobe filling helium donor, an accretor mass 
Mi = 1.3 M would be required. This is incompati- 
ble with our mass ratio measurement. Note that mass 
transfer reduces the orbital period change rate relative 
to the case of a detached binary, so that higher masses 
would be required; but in practice, it is likely that the 
mass transfer is not yet sufficiently developed to have 



a big influence (Deloye et al. 2007). We can conclude 
that the donor star is probably at least twice as mas- 
sive as a fully-degenerate helium donor would be. For 
masses M 2 = 0.27 M , Mi — 0.55 M Q and assuming 
that the radial-velocity semi-amplitude of the wobbling 
Hen 4686 line represents the projected orbital velocity 
of the accretor, we obtain an inclination i sa 38° of the 
orbital plane. 

4.3. HM Cnc as a Gravitational Wave Source 

With the system parameters derived in the previous 
section, we can calculate the gravitational-wave strain 
amplitude at Earth. The distance to HM Cnc is the 
largest remaining uncertainty; it is probably ~5 kpc 
based on the expected sizes an d measured temper atures 
of the white dwarf components ( Barros et al.|2007 |. This 
gives a dimensionless gravitational- wave strain ampli- 
tude h ~ 1.0 x 10~ 22 . This, together with the design sen- 
sitivities of space-borne gravitational wave detectors like 
the Laser Interferometer Space Antenna (LISA), makes 
HM Cnc one of the easie st detectable sou rces of gravita- 
tional waves known (see Nelemans||2009 ). 



We thank Susana Barros for helping us with 
the ephemeris of HM Cnc. GHAR was sup- 
ported by NWO Rubicon grant 680.50.0610 to G.H.A. 
Roelofs. AR acknowledges support through NASA grant 
NNX08AK66G. Based on data obtained at the W. M. 
Keck Observatory, which is operated as a scientific part- 
nership among the California Institute of Technology, the 
University of California and the National Aeronautics 
and Space Administration. The Observatory was made 
possible by the generous financial support of the W.M. 



REFERENCES 



Barros, S. C. C, Marsh, T. R., Groot, P., Nelemans, G., Ramsay, 
G., Roelofs, G., Steeghs, D., Wilms, J. 2005, MNRAS, 357, 
1306 

Barros, S. C. C. et al. 2007, MNRAS, 374, 1334 

Cropper, M., Harrop-Allin, M. K., Mason, K. O., Mittaz, J. P. 

D., Potter, S. B., Ramsay, G. 1998, MNRAS, 293, L57 
Dall'Osso, S., Israel, G.L., & Stella, L. 2006, A&A, 447, 785 
Dall'Osso, S., Israel, G.L., & Stella, L. 2007, A&A, 464, 417 
D'Antona, F., Ventura, P., Burderi, L, & Teodorescu, A. 2006, 

ApJ, 653, 1429 
Deloye, C. J., & Taam, R. E. 2006, ApJ, 649, L99 
Deloye, C. J., Taam, R. E., Winisdoerffer, C, Chabrier, G. 2007, 

MNRAS, 381, 525 
Hakala, P., Ramsay, G., Wu, K., Hjalmarsdotter, L., Jarvinen, S., 

Jarvinen, A., Cropper, M. 2003, MNRAS, 343, L10 
Hakala, P., Ramsay, G., & Byckling, K. 2004, MNRAS, 353, 453 
Israel, G.L. et al. 2002, A&A, 386, L13 

Israel, G.L. et al. 2004, Mem. Soc. Astron. Ital. Suppl., 5, 148 
Marsh, T. R., & Home, K. 1988, MNRAS, 235, 269 
Marsh, T. R., & Steeghs, D. 2002, MNRAS, 331, L7 
Munoz-Darias, T., Casares, J., & Martmez-Pais, I.G. 2005, ApJ, 
635, 502 

Motch, C, Haberl, F., Guillout, P., Pakull, M., Reinsch, K., & 

Krautter, J. 1996, A&A, 307, 459 
Nelemans, G., Portegies Zwart, S.F., Verbunt, F., & Yungelson, 

L.R. 2001, A&A, 368, 939 
Nelemans, G. 2009, Class. Quantum Gravity, 26, 094030 
Norton, A.J., Haswell, C.A., & Wynn, G.A. 2004, A&A, 419, 1025 



Oke, J.B. et al. 1995, PASP, 107, 375 

Ramsay, G., Cropper, M., Wu, K., Mason, K.O., Hakala, P. 2000, 

MNRAS, 311, 75 
Ramsay, G., Hakala, P., & Cropper, M. 2002a, MNRAS, 332, L7 
Ramsay, G., Wu, K., Cropper, M., Schmidt, G., Sekiguchi, K., 

Iwamuro, F., Maihara, T. 2002b, MNRAS, 333, 575 
Reinsch, K., Steiper, J., & Dreizler, S. 2007, in ASP Conf. Ser. 

372, Spectral Analysis and Constraints on the Nature of the 

Ultra-Compact Binary RX J0806. 3+1527, ed. R. Napiwotzki & 

M.R. Burleigh (San Francisco, CA: ASP), 419 
Roelofs, G.H.A. et al. 2009, MNRAS, 394, 367 
Steeghs, D., Marsh, T. R., Barros, S. C. C, Nelemans, G., Groot, 

P. J., Roelofs, G. H. A., Ramsay, G., Cropper, M. 2006, ApJ, 

649, 382 

Strohmayer, T.E. 2002, ApJ, 581, 577 
Strohmayer, T.E. 2003, ApJ, 593, L39 
Strohmayer, T.E. 2004, ApJ, 610, 416 
Strohmayer, T.E. 2005, ApJ, 627, 920 
Strohmayer, T.E. 2008, ApJ, 679, L109 

van der Sluys, M. V., Verbunt, F., & Pols, O. R. 2005, A&A, 440, 
973 

Wood, M.A. 2009, MNRAS, 395, 378 

Wu, K., Cropper, M., Ramsay, G., Sekiguchi, K. 2002, MNRAS, 
331, 221 

Yungelson, L.R. 2008, Astron. Lett., 34, 620 



